Stripes, pseudogaps, and SO (6) in the cuprate superconductors 
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We briefly summarize two related calculations. First, we 
demonstrate that the instabilities (either nesting or pairing) 
associated with the high-Tc cuprates can be described by 
an 5*0(6) transformation group. There are two independent 
6-dimensional representations ('superspins'). One superspin 
combines Zhang's 5-component superspin with a flux phase 
instability; the other involves a charge density wave, s-wave 
superconductivity, and an exotic spin current. 

The second calculation is a self-consistent slave boson cal- 
culation, which provides a good description of the doping 
dependence of the photoemission dispersion in terms of dy- 
namic striped phases. The stripes are stabilized by strong 
electron-phonon coupling, and provide evidence for a doping- 
dependent crossover between the two superspin groundstates. 



I. SO{<o) INSTABILITY GROUP 



+1 



+1 - 



- 



-1 - 





4 




OSDW- 


OsDWz- 
-OjC 


OsDW+ 




Ad 





+1 









Ojs- 




Ojs+ 




As 





> 

m 
T— I 

o 



I 

o 
o 



X 
S3 



ir]- 








Ilx+ 










-Ay 








-A, 








Q 


V+ 


-aix- 


-iUy- 





FIG. 1. Matrix Representation of S0(6), using the short- 
hand 0± = O ± 

Zhang |l| has proposed that the physics of doping the 
cuprate superconductors is controlled by an 50(5) trans- 
formation group, having a 5-dimensional superspin rep- 
resentation which mixes spin-density wave (SDW) and 
d-wave superconducting components. This is in contra- 
diction to earhcr findings that in the Hubbard model 
at half filling d-wave superconductivity is symmetry- 
equivalent to a fiux phase Here, we show that Zhang's 
50(5) group is a subgroup of an 5*0(6) group, the in- 
stability group of the two-dimensional (2D) Van Hove 
singularity (VHS) The corresponding 6-component 
superspin combines Zhang's superspin with a fiux-phase 
instability. 

Consider an electronic state with a two-fold orbital de- 
generacy, ipi:''p2- Including spin and charge conjugation. 
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FIG. 2. S'0(4) weight diagram of (top) and Vn (bot- 
tom), where Sz (Jz) is the z component of spin (pseudospin) . 

8 single particle operators can be formed, or 8 x 7/2 = 28 
pair operators (taking account of the Pauli exclusion 
principle). These operators can be combined into a 
representation of an 50(8) group, the pair group jj] 
which controls the various pairing and nesting instabili- 
ties which can arise. For the quasi-2D cuprates, the or- 
bital degeneracy is associated with the VHS's along (tt, 0) 
and (0, tt); because of the large density of states (dos) at 
a VHS, the Fermi surface can be approximated by just 
two points [^y^, similar to a one-dimensional metal. 

Consider a supercell of size 2a x 2a, each cell containing 
four atoms. Let j be a lattice vector of the original lattice, 
and r be a lattice vector of the supercell lattice. If the 
Cu-site operators are a's, then the VHS operators are 



il'2n('r) — —(a^ + a^, - — a^, ~ — a^, ~ , . ). (I) 

2 i+x,{T i+y.cr i-\-x-\-y,{T ^ V / 

(The remaining operators are located in other parts of 
the Brillouin zone, and play a negligible role in nesting 
phenomena.) 

A special role is played by the operator t = 
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J2'Kai'^L(.r)i'ia{r) - ^l„{r)^p2cr{r)), which sphts the 
VHS degeneracy. The operators which commute with 
T form an 50(6) transformation group, Fig. 1, with two 
fundamental 6-dimensional representations, Fig. 2. The 
momentum-space representation of the superspin opera- 
tors is shown in Table I; the full list is in Ref. Q. Since 
the generators of SO (6) can be represented as antisym- 
metric 6x6 matrices. Fig. 1 assigns each operator to the 
position of the corresponding non-zero matrix element, 
to produce the correct SO{6) Lie algebra 



(2) 



Note that by deleting the second row and column, an 
SO{5) subalgebra is generated, which is equivalent to 
the one postulated by Zhang [Q. 

In Fig. 2, the pseudospin operator m is 



J- = 277, J+ = 2?7^ J. = Q, 



(3) 



and OsDW± = TiOsDWx ± iOsDWy)/V^- The su- 
perspin Vg consists of a spin-density wave spin triplet 
(OsDw) and a pseudospin triplet combining d-wave su- 
perconductivity {Ad, A^) and an orbital antiferromagnet 
(Ojc) equivalent to the flux phase operator Q. The su- 
perspin Vh consists of a spin current spin triplet (Ojs) 
and a charge-density wave (Ocdw), s-wave superconduc- 
tivity (As, A|) pseudospin triplet. These operators are 
equivalent to those introduced by Schulz Q and, in a 
related context, by Halperin and Rice 

Table I: Fourier Transforms of Superspins 
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In weak coupling, the instabilities are controlled by 
their respective susceptibilities and the form of the quar- 
tic interaction terms. For an energy dispersion 

— — 2t(cos kr^a + cos kya) + At' cos k^a cos kyU, (4) 

then when t' = and at exactly half filling (square Fermi 
surface), the model has an extra pseudospin symmetry 
[0, and all the bare susceptibilities are equal Q. The 
dressed susceptibilities are 



Xoi 



where Gi is related to the quartic coupling. Since all 
12 Xoi's are equal (at half filling), the ground state in- 
stability is controlled by the interactions Gi. For a pure 
Hubbard model, the U term breaks the 50(6) symmetry, 
leading to a negative Gsdw and a positive Gcdw, Gs, 
Fig. 3. Thus at half filling the ground state is antiferro- 
magnetic, Osdw- 
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1 + 2xo^G, 



(5) 



FIG. 3. Splitting of 50(6) superspin degeneracies by inter- 
action U and by doping x. 

Doping away from half filling breaks the pseudospin 
symmetry, leading to two distinct susceptibilities, xoo 
for the particle-hole operators and xo2 for the particle- 
particle operators, Fig. 4. Since the susceptibilities are 
enhanced by nesting, their magnitudes fall off with dop- 
ing. However, the magnitude of xoo falls much more 
rapidly, leading to a possible scenario for d-wave super- 
conductivity. Fig. 3 IJ]: if there is an attractive interac- 
tion for d-wave superconductivity, it will be overwhelmed 
by the strong divergence in xs dw at half filling, but with 
doping, the SDW instability falls off much faster, lead- 
ing to a possible crossover to d-wave superconductivity. 
This is consistent with earlier renormalization group re- 
suhs 

In a more realistic, strong coupling model, a number of 
additional features must be considered. First, additional 
hopping parameters, such as the second-neighbor t' term 
in Eq. ^, also break the pseudospin degeneracy, Fig. 5. 
Note that the susceptibility diverges away from half fill- 
ing, at the VHS, but even at optimal doping xo2 >> Xoo- 
Moreover, by comparing Fig. 4 it can be seen that X02 
is actually larger, at optimal doping, for t' ^ Q. Since 
strong coupling pins the Fermi level near the VHS over 
an extended doping range ||ll|,|2|, this enhancement of 
Xo2 can greatly enhance the possibility for d-wave super- 
conductivity to arise on doping. 

There is considerable evidence for phonon coupling and 
structural instabilities in the doped material pO| ]. Hence, 
it is interesting to see whether strong electron-phonon 
coupling can cause a crossover to a groundstatc involv- 
ing V/j. Since the CDW operator is proportional to the 
difference in hole population of the odd and even sublat- 
tices, it must vanish in the strong coupling limit at 
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FIG. 4. Susceptibilities xoo (dotdashed line) and xo2 (solid 
line) and density of states (dashed line) vs. band filling 1 + x 
for the one band model with t' = 0. 




FIG. 5. As in Fig. g but with 2t' = -0.3t. 

half filling, when there is one hole per Cu. For finite dop- 
ing there can be a CDW order, with < Ocdw >oc x, if 
all the holes are confined to one sublattice. The resulting 
structure is equivalent to that found in La2-2:Sr2:Ni04 
for X - 0.5 0]. Moreover, such a structure is closely 
related to the holon condensate introduced to over- 
come parity-violation problems in anyonic models of su- 
perconductivity. In these models, two species of holons 
are introduced, of opposite parity violation. A coherent 
superposition of the two preserves parity, while produc- 
ing a holon condensation, wherein all holons are confined 
to one sublattice, leading to a real-space CDW [Tsl] . 

In the following section, we describe a slave boson cal- 
culation jl^ which self-consistently takes into account 
both strong coupling effects and strong electron-phonon 
coupling. We show that VHS pinning can provide a nat- 
ural explanation for the observed stripe phases [^jl7j 
in these materials, as well as pseudogaps and extended 
VHS's. 



II. CDW ORDER AND PHASE SEPARATION 

In the strong coupling limit, the dominant physics is 
quite different. At exactly half filling, there is a metal- 
insulator transition that is unrelated to nesting: a Mott 
transition, or charge transfer insulator transition in the 
three-band model of the cuprates which will be employed 
here. The nesting/pairing instabilities arise on top of 
this transition, and are strongly constrained by the large 



Hubbard U, which limits double occupancy of the Cu 
orbitals. In the slave boson approach employed here, the 
limit J7 — > oo is taken (no double occupancy). To include 
spin effects, it is necessary to go beyond the standard 
slave boson approach, either by introducing four slave 
bosons iQ or by introducing a Heisenberg exchange J 
I p^ - the latter approach is followed here. 

Thus, the starting point of the analysis is the three- 
band model Hamiltonian: 

^ = E(^4,a«'j^- + E tcuo[dl,P,+i^ + (c.c.)] 

S' 

+jT.<^d,^^"^+..d^+,J. (6) 

In this equation, (p^) is a creation operator for holes 
on Cu (O), j is summed over lattice sites, 6 over Cu-0 
nearest neighbors, 6' {S") over 0-0 (Cu-Cu) neighbors, 
and c.c. stands for complex conjugate. The band param- 
eters are A, the splitting between the Cu and O energy 
levels, tcuO, the Cu-0 hopping parameter, too, the 0-0 
hopping parameter, and J the exchange constant. 

In the above equation, the limit U oo has been 
taken, and a slave boson introduced to prevent Cu double 
occupancy. The parameters tcuO ^^nd A are renormalized 
values, to be solved self-consistently jlj]. In the exchange 
energy, the d'^-term is decoupled in mean-field as 



j,6" .a.a' 



J 



i,<5",o" 

where Ns is the number of unit cells and 
A„ = < d,^d]^ >= Aie'^' 



(7) 



(8) 



Depending on the choice of phase, a variety of magnetic 
states are possible. Here the only choices considered are 
6ij = 0, corresponding to a uniform phase (only short- 
range magnetic correlations), and 9ij — ±7r/4 for a flux 
1^ phase, with the ± sign chosen so that the net phase 
change around any plaquette is tt. 

Finally, electron-phonon coupling is included by the 
substitution tcuO = ^(1 i S), with S due to a CDW dis- 
tortion: for a breathing mode, all four Cu-0 bonds are 
equivalent for any given Cu, but all are long (tcuO — 
t{l — 5)) for one sublattice, and short for the other. Note 
that this model includes only two of the twelve possible 
superspin components, but already gives rise to a com- 
plicated phase diagram. It should describe the Vg Vh 
(Ojc Ocdw) crossover. 
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A charge transfer insulator transition is found at half 
filling when A exceeds a critical value . Near half fill- 
ing, the flux phase is stable, and the CDW is suppressed 
by strong coupling effects. However, if the electron- 
phonon coupling is strong enough, there is a crossover in 
the doped material to a CDW phase, which is maximally 
stable near the doping at which the Fermi level crosses 
the bare VHS. As had been previously noted |Q, the 
flux phase provides a good description of the measured 
energy dispersion in insulating cuprates |22|, Fig. ^, while 
the CDW fits the dispersion near optimum doping [ p3[ , 
in particular reproducing the extended VHS. [Note how- 
ever that the present calculations underestimate the total 
dispersion by a factor of 2.] 
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FIG. 6. (a) Energy dispersion for the flux phase at half 
filling Xi — 0^ (dashed line) and the paramagnetic phase at 
X = 0.22, Vep = 0.6eV (solid line), close to the minimum free 
energy, and for the fluctuating stripe phase model, for Uc — 
0.5 (dotdashed fine) and 0.75 (dotted fine). Data from under- 
doped Bi-2212 (diamonds and x's) ||| or SCOC (squares) 
[ p2[ are plotted as E/2. Special points of the Brillouin zone 
are X — (7r,0), S = {n,n). (b) Density of states vs. energy, 
for the same parameters. 

The situation for intermediate dopings is more com- 
plex. The free energy minimum at finite doping associ- 
ated with the CDW leads to an instability toward phase 
separation between magnetic (flux phase) and charged 
(CDW) domains. There is considerable experimental evi- 
dence for such phase separation in the cuprates ||2j,^,|l^ , 
which can be macroscopic if the doping counterions are 
mobile (e.g., the interstitial oxygen in La2Cu04+5 [p6| . 
More commonly. Coulomb repulsion between holes re- 
stricts the phase separation to a nanoscopic scale [ p7| , 
typically forming striped phases |l^ , p^ . The driving 
mechanism for this phase separation has variously been 
ascribed to large J [^ , p9| , to the nearest neighbor 
Coulomb energy V |3C|], or, as in the present case, to 
VHS-induced CDW formation Q. 

In modeling the striped phase, a static stripe pattern 
produces too much structure in the dispersion, due to 



minigaps . If a fluctuating striped phase is approx- 

imated by taking the weighted average of the end-phase 
band parameters, good agreement with the experimen- 
tal dispersion is found. Fig. ^ In this figure Vc = x/xc 
is the fraction of charged stripes, with Xc the doping of 
the uniformly doped charged end phase, here taken as 
Xc — 0.288. The same model also describes the crossover 
from small (hole pockets) to large Fermi surface with 
doping. Fig. 
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FIG. 7. Fermi surfaces for the fluctuating stripe phase 
model, for i^c = 1.0 (dotted fine), 0.5486 (solid line), 0.5 
(dashed line), and 0.25 (dotdashed line). 

A number of points should be made. (1) The split- 
ting of the bands near the X-point, (tt, 0), provides a 
good descrpition of pseudogap formation in the cuprates 
- not only the photoemission data. Fig. ^, but also ther- 
modynamic data from heat capacity and magnetization 
1 32 1, as well as the temperature-doping phase diagram 
|33|. The smooth doping dependence masks a crossover 
from spin- like to charge- like instability - from Vg — > Vh. 
Accordingly, the underdoped results near half filling are 
in good agreement with recent Monte Carlo calculations 
Q. (2) From Fig. it is clear that the pseudogap 
is directly associated with the splitting of the VHS de- 
generacy, and hence with the properties of 50(6). (3) 
Campuzano | |35[ | reports a different doping dependence of 
the (tt, 0) photoemission, namely that there are two inde- 
pendent features, one a sharp quasiparticle peak which 
is near the Fermi level at optimal doping, and broad- 
ens severely in underdoped samples, and the second a 
broad peak which is already present at ~200meV below 
the Fermi level in optimally doped material, and grad- 
ually shifts to 300meV with increased underdoping. A 
dynamic fluctuation model can accomodate this in terms 
of a dynamic average of the two end phases. (4) There is 
considerable evidence for short-range CDW order in the 
doped cuprates, summarized in Ref. jl^. Section 9.2, of- 
ten associated with CuO bond stretching p&. (5) In the 
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present model, the pseudogap is associated with mag- 
netic or structural instabilities which compete with su- 
perconductivity. Hence, there should be a separate su- 
perconducting gap, which scales with and This 
gap may be related to the spin gap and resonance peak 
seen in neutron scattering measurements of the magnetic 
susceptibility near S = (tTjTt) (6) In the overdoped 

regime, the CDW undergoes a quantum phase transition, 
TcDW 0. It has been suggested that such a QCP can 
lead to high-Tc superconductivity [|8|. 
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